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We study meson-meson scattering in a soluble model which describes asymptotically 
free mesons and confined quark-antiquark pairs via coupled channels. Concretely, the two 
scattered mesons are assumed to interact through s-channel meson-exchange diagrams. Fur- 
thermore, we develop a perturbative expansion of the model, and show that the thus found 
resonance pole positions, including contributions up to fourth order in perturbation theory, 
completely fail to reproduce the exact results. This shows that the resonance predictions 
| based on perturbative approximations in quark models may be highly unreliable. 
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§1. Introduction 



Hadron spectroscopy in present times appears to have become subdivided into 
many different fields of interest, like conventional quark spectroscopy, hadronic 
molecular states, dynamically generated resonances, tetraquarks and pentaquarks, 
glueballs, gluonic hybrids, and so forth. The advancement of detector and analy- 
sis techniques at the many new experimental facilities in the intermediate-energy 
range has been resulting in the observation of more and more hadronic states, many 
of which no not seem to fit into the traditional quark spectrum of qq mesons and 
' qqq baryons. This has led to the investigation of other possible configurations that 

might be viable within QCD. On the one hand, this makes it important and inter- 
esting to explore the structure and properties of hadrons, which may shed light on 
the dynamics of strong interactions at low energies. On the other hand, though, 
a very careful analysis is required, in order to avoid confusions and controversies. 
We are convinced that the most important need of present-day research in hadron 
spectroscopy is some consensus on the ideal quark model to confront with the data. 
And right now is probably the best moment to tackle this problem, since quite a lot 
of data, at least concerning meson spectroscopy, is being produced (see for example 
RefsPJlSffl©.©©,©®^ and even more data, of better quality, is expected in the 
near future. All this could help gather a vastly improved understanding of quarkonia 
and other mesonic resonances. However, in such a situation it is extremely impor- 
tant to first agree on the hadron spectrum, as obtained from a chosen quark model, 
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since only then valid conclusions can be drawn about possible incompatibilites with 
standard quark configurations. 

Indeed a lot of work is being done on trying to understand the conventional 
and unconventional structure of different "exotic" hadrons. For example, a very 
interesting analysis of the quark content and possible molecular nature of many 
newly found mesonic states has been carried out on the basis of QCD sum rulesp^'HJ 1 
while an attempt to distinguish between quarkonium and hybrid states was made in 
RefPJ Another very interesting investigation, namely of the quark and molecular 
content of baryon resonances, was reported in RefP^ 1 Furthermore, RefP nicely 
explained the phenomenon of dynamically generated resonances and the concept of 
dynamical reconstruction, in order to be able to pinpoint (dominantly) qq states. 
All these meticulous efforts and the corresponding results would get more merit, if 
a universally agreed conventional hadron spectrum was known. 

The main issue is that experimental evidence for possible resonances is obtained 
from total or partial-wave cross sections, as well as angular distributions and de- 
cay modes. In order to interpret the data, one needs a model, since perturbative 
QCD cannot be used at low energies. In the present paper, we will focus on non- 
exotic meson-meson scattering. Nevertheless, the results can easily be generalized 
for application to hadron-pair production.^ We intend to describe the cross sec- 
tions and resonance pole positions for meson-meson scattering in an as large as 
possible energy range, rather than focusing on just one peak. At this point enters 
the main philosophy behind the model, namely that the enhancement structure of 
cross sections in non-exotic meson-meson scattering stems from the quark-antiquark 
spectrum. Consequently, for such reactions we must study a coupled system of a 
qq state and non-exotic two-meson channels. This is an absolutely minimal require- 
ment for modeling the cross sections in this case. Further extensions to multiquark 
or hybrid resonances might be contemplated in case the minimal model turns out 
not to reproduce the experimental data sufficiently well. Furthermore, the proposed 
strategj^D 1 '^ simultaneously covers two-meson molecules and qq systems. Hence, 
the physical solutions are not just pure qq or molecular states, but rather mixtures 
of these two configurations. One could then try to find out which component is 
dominant QHJiEOJSSJ However, that lies beyond the scope of this paper. In more 
popular terms, one might refer to the states obtained from such a model as qq sys- 
tems surrounded by a meson cloud. In the past, our approach was also called the 
unitarization scheme for the quark-antiquark systemPU'EZJ'EBJSJ 

The model we are going to use treats confined quarks and hadronic decay 
channels on an equal footing, via coupled channels, regardless of whether the en- 
ergy is above or below the thresholds of the decay channels. It was developed 
in Kefs-p^'^'^'^'^ and has been extensively used to study the properties of 
mesonic resonances (for some of the recent works, see Ref s PSJ 'EHJ 5 l3oJ ,[3T]i ,[32j ,[33]i anc j 
references therein). The effective meson- meson potential in the model consists of 
s-channel exchange of a confined qq pair, with radial quantum number running from 
to infinity and orbital angular momentum compatible with total J PC . The im- 
portance of hadron dynamics in understanding the meson spectrum for low and 
intermediate energies already becomes obvious from the fact that it easily gets en- 
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ergetically favorable for the "string" between the quark and the antiquark to break, 
alongside the creation of a new and light qq pair, which then may lead to hadronic 
decay. An even stronger indication comes from the light scalar mesons, whose un- 
conventional nature is related to their very strong coupling to S-wave two-meson 
channels P^l The latter work also shows that the inclusion of bare- meson exchange 
in the s-channel, in addition to the meson-meson contact interactionJ^'^'^'^'ED 1 
in certain cases leads to finding new states which might be absent when considering 
contact interactions only. This phenomenon was explained very neatly in Ref.p3 
where it was shown that, for the study of a restricted energy range corresponding 
to a particular resonance, the contribution from different diagrams involving meson 
exchange with different quantum numbers gives rise to a constant interaction, which 
is equivalent to considering a contact interaction in unitarized models. It was further 
shown in RefP3 that, in order to understand a larger energy range, covering several 
resonances, meson-exchange diagrams are required as well. This explains why the 
common use of contact interactions in unitarized models, to study dynamically gen- 
erated hadron resonances, works quite we l\^MBEBE2bE!BEQEBMbM$ However, 
from Refp3 it becomes clear that the development of a broader and more general 
perspective for hadron spectroscopy requires the treatment of quarks and hadrons as 
coupled systems. In this paper, we will show that not only the handling of coupled 
mesons and quarks is necessary, but also the full solution of the scattering equations 
is essential. In particular, we will demonstrate that approximating resonance pole 
positions perturbatively leads to unreliable results. 

In the next section, we will first describe the exact formalism, followed by the 
construction of a perturbative expansion thereof. In the subsequent sections, we will 
choose some specific examples to show that no meaningful results can be obtained 
from the perturbative expansion. Finally, we will summarize the detailed discussions 
in the paper. 

§2. Formalism 

We study meson-meson scattering using a model in which quarks and hadrons 
are considered coupled systems. The formalism amounts to solving a scattering 
equation for mesons, with the lowest-order term of the Born series given by an 
effective interaction due to the exchange of a confined qq pair. The potential between 
the latter pair is written in terms of a harmonic oscillator. The eigenenergies of the 
harmonic oscillator thus correspond to the bare qq spectrum. The effective meson- 
meson potential, which is the lowest-order term of the full scattering amplitude, 
involves an infinite sum over all the diagrams with s-channel qq exchange, having 
different radial quantum number n (0 < n < oo), as shown in Fig. [1] below. 

Although not strictly necessary, it is illustrative to consider a formulation of 
the model in terms of a coupled system of nonrelativistic Hamiltonians. However, a 
rigorous derivation in terms of a sum of meson loops, which leads to the same final 
result, is also possible. 

Since the model is based on coupling a confined quark-antiquark pair and a 
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Fig. 1. The meson-meson potential, i.e., the Born term of the full scattering amplitude, which 
involves the exchange of a confined qq pair, with radial quantum number n running from to 
oo. 

meson pair, we describe the system by the equations 

ffc^c(r) + V T (r)iP f (r) = Eip c (r), (2-1) 
Hf^f{r) + V T (r)Mr) = E^ f (r), (2-2) 

where the subscripts "c" and "/" (here and throughout this article) refer to the 
confined quarks and free mesons (i.e., considering them plane waves), respectively, 
and Vt is the transition potential between the two sectors. 

H c and Hf describe the Hamiltonians of these sectors, reading 

H c = -^+m q + mg + V c (r), (2-3) 

H f = -^ + M 1 +M 2 , (2-4) 
where the confining potential is assumed to be a harmonic oscillator, viz. 

V c = ^ c co 2 r 2 , (2-5) 

with fi c and ui the reduced mass and frequency of the qq system, respectively. Fur- 
thermore, Mi, M2, and m q , mq are the meson and quark masses, respectively. Our 
choice of a harmonic-oscillator potential in the confined sector is based on earlier 
observations of regular spacings in the quarkonium spectrap^'^ which seem to 
be confirmed by states found in even the most recent experiments 'SHJ ESJEB j n 
writing down the equations above, we have assumed only one qq and one meson- 
meson channel, for the sake of simplicity. These equations can be straighforwardly 
generalized to the multichannel case, taking then a matrix formPD 
Now, Eqs. (|2-2p and (|2-ip can be rewritten as 



(E - H c )Mr) = V T (f)^ f {r), (2-6) 
(E-H f )i; f (r) = V T (r)Mrl- 

Then, the confinement wave function ipd'r) must be eliminated from the equations, 
as it never develops into an asymptotic state. Thus we get 



V7(f) = {E- H f )~ l V T {E - H c )- l V T ^ f {r). 



(2-7) 
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From first principles of standard scattering theory, we can conclude that the factor 

V T (E - H C )- X V T (2-8) 
acts like an "effective" meson-meson potential, which, if denoted by Vmm, implies 

(Pf | Vmm I P' f ) = (P f | V T (E{Pf) - H^Vt \ P' f ), (2-9) 
where the total center-of-mass energy (CM) E is given by 

E(p f ) = ^^+M 1 + M 2 , (2-10) 

with fif the reduced mass of the two mesons, and Pf (P'f) denoting the CM momen- 
tum of the two-meson initial (final) state. 

Furthermore, we denote the energy eigenvalue of H c by E n i, i.e., 

E nl = uj{n c + l c + 3/2) + M q + M q , (2-11) 



and the corresponding eigensolutions by (f c \ n c , l c , m c ). By introducing in Eq. (|2-9p 
a complete set corresponding to this state, we get 

{Pf I Vmm I P'f) 

= Y, (Pf\VT\n c ,l c ,m c )(n c ,l c ,m c \(E(P f )-H c )- l V T \P'f) 

= E (Pf\VT l ^^^ cl V T \P' f ), (2-12) 

which, upon further introduction of several complete sets corresponding to the 
meson-meson configuration space, gives 



(P f | Vmm I P'f) = £ J d\f J afr'f j d\" f j 



n c ,l c ,m c " (E(Pf)-Enl) 

x(f f | V T I r))(r] | n c ,l c ,m c )(n c ,l c ,m c \ rj'){rj' | V T \ | P' f ). (2-13) 

For the transition potential, we take a local delta-shell function of the form 

(f f | V T I f* f ) = S(r f - a)5 3 (r f - r'f). (2-14) 

This form of potential has been proven useful in describing the breaking of the 
color stringJED The A and a in Eq. (|2-14p are the two parameters of the model, with 
the former being the coupling of the meson channel to the quark channel, and the 
latter an average distance between the quarks. The coupling A is varied between 
and 1 in the present study, with A = corresponding to decoupled meson and quark 
systems. Since the meson-meson state is considered a plane wave, decoupling would 
result in a pure ( "bare" ) qq spectrum. On the other hand, A > 1 represents strong 
coupling to the meson-meson channel. The parameter a is taken in the range 3-5 fm. 
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Using the above form for Vr, and the normalization (r | p) = e % ^' r /(2tt) 3 / 2 , 
Eq. (|2-13p becomes 



(Pf I V M m I P' s ) = £ / -#5 I4^*- P > f >^t 



a x 



(ff I n c ,L,m c )(n c ,l c ,m c I r^) X . 
x W ' ^ A - f —5(r' f - a)e lP f r f. (2-15) 

- ^ Vca 1 

The functions (r} | n c ,l c ,m c ) = Ym f (f)(\ff \ \ n Cl l c ,m c ) represent an overlap wave 
function of the meson-meson — > qq vertex. In order to determine this overlap func- 
tion, we assume a mechanism for the transition from the meson channel to the 
quark channel, or vice versa. To describe this mechanism, let us consider a confined 
quark-antiquark pair with total spin, angular momentum, intrinsic spin, and radial 
quantum number j, l\, si, and m, respectively. As one usually expects from QCD 
at low Q 2 , it is assumed that the string between the quark and the antiquark breaks 
with the creation of a new qq pair. This new pair is assumed to have vaccuum 
quantum number, corresponding to a 2s2+l hj 2 = 3 Po state. This yields two qq pairs, 
which then rearrange to produce two mesons with quantum numbers j[ , l[ , S[ , n\ 
and j' 2 , 1' 2 , S' 2 ,n 2 , respectively (see Fig. [2]). 



qq M_ „ j'i> !'r s \> n \ 

j i555555*' Ti = 



M ~i'r V t s 'r n ' 2 



} 




M 

Fig. 2. The transition vertex from the quark channel to the two-meson channel. 

Mathematically, the rearrangement of the two qq pairs, one of which is a 3 Po 
state and the other having the quantum numbers of the decaying meson, can be 
expressed by treating them as four independent nonrelativistic harmonic oscillators. 
Let us label them with numbers, and assume that system (1+2) represents the quark- 
antiquark pair under consideration, while system (3 + 4) stands for the newly created 
3 Po q<l pair. This can be treated as a four-body problem, which can be reduced 
to a three-body problem in the global CM system, by considering the coordinates 
(momenta) of the CM of the (1 + 2), r\ 2 {p\ 2 ), and (3 + 4), ^34(^34) systems, along 
with their relative motion ^1234(^1234)- The situation after the transition is described 
by assuming that the (1 + 4) and (3 + 2) systems represent the two meson state. 
This can be represented diagrammatically, as shown in Fig. [3l 

The Hamiltonians for the systems, before and after the transition, can be written 
in the global CM frame as 

^before = 7^{ r 12 + r 34 + r 1234 + Pl2 + ^34 +P1234}, (2-16) 
Rafter = 7^{ r 14 + r 32 + r 1432 + Pl4 + #32 + P1432}, (2'17) 
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Fig. 3. The four quarks, labeled by numbers, before and after the transition. It is assumed that 
before the transition the string between a qq pair (1 and 2) breaks to form a new pair (3 and 
4), with 3 Po quantum numbers, and the whole system rearranges so as to produce two mesons 
made of the quark pairs 1-4 and 3-2. 



where r^- (pij) is the coordinate (momentum) of the ij system, and r^jy (pijkl) is the 
relative coordinate (momentum) of the CM of the ij and kl systems. The transition 
of e.g. the coordinates of the four quarks can be expressed in terms of an orthogonal 
transformation matrix a (as explained in Ref.^J), i.e., 



/ >| 




-32 




\ H432 J 





"11 "12 "13 
«21 "22 "23 
"31 "32 "33 




(2-18) 



Thus, if the wave functions of the systems before and after the transition are 
written as V , f l) / )m }('"i2, -34, 7*1234) and xf n /^ iTO /}(n4, -32, n.432), respectively, then 



(Vf„,i, m} (n2,r 3 4,n234)| = J dr i4jdr 23 J 

n'.l'.m' 



dr 



1234 



(^{n,l,m}( ri2 ' rU ' n234 ) I *{n',l',m'}( r l4, r 32 , r M32 )) {X{ n ' ,V ,m>} (-14, -32,-1432) 



Now we define a transformation matrix Df; ^ 



{{n,l,m},{n',l',m'}} 



HS 



D 



{{n,l,m},{n',l',m'}} 



dr u J dr 23 J dri 2M (^f n ^ m} (r 12 ,r u ,ri 2U ) \ X{ n >,l',m'}( r i4> -32, -1432)), 
for which the following analytical expression was obtained in Ref.P^ 



V 



{{n,t,m} , { 
4 



n'J'm'}} (W ' I 1 "'}) 



'}} 

N 

IJ(-1)»* 



(2-19) 
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TV TV 

nn 

i=i j=i 



Uij + 1 



r (mj + i)r (n i:j + Uj + 1) 



(2-20) 



We denote the elements of this transition D-matrix by g n , which give the proba- 
bility of the transition from a particular meson channel to a particular quark chan- 
nel, and vice versa. These g n precisely correspond to the overlap wave functions 
yfaj fi c (\ff \ | n c , l c , m c ) in Eq. (|2-15|) . Replacing these overlap functions by n c 9n/ (v^) 
turns Eq. f|2- 15[) into 



(Pf | V MM I P' f ) 



/ df2 rf / ain rf e- lP f af f 



A \ Vc9n 



ji c a J a 



Y^ f (f f )Y^ f (ff') ip,.^. 



E(P f ) - E nl 



\ 2 a 
(2vr) 3 



£ J dQ rf J df2 r , f e- iP r af fY% f (f f )YX (r/0 x 



9n 



E(P f ) - E nl 



iP't-ar', 
-e f 1 . 



Using the standard relations for spherical harmonics 

j df2 rf e- iP ^f a Y^(r f ) = {-if^j^P^Y^Pf) 
in Eq. (l2 T 2T|) . we get 



(2-21) 



(2-22) 



{Pf I V MM I P'f 
\ 2 a 



(2vr) 3 
(2vr) 3 



^ (-i)'/47riK^/a)^ / / (P / )(i)^47rj / ( J P}a)y^/( J P}) 



9l 



£ H 2 )" {^?3l (Pfa)ji {P' f a)Yk f {Pj)Y*j (p' f ) 



E(P f ) - E n i 
9l 



E(P f )-E nl 
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X 2 a 
(2vr) 3 



^(-i 2 ) ; /(4vr) 2 iKP / a)iK^a) : 



2l f + l 



7l c .l c 



4tt 



9n 



E{P f ) - E nl 



2vr 2 



J2(2i f + i)F lf (p f . pi) 3l{Pfa ) jl{ pi a ) j2 



lc = 



n c =0 



E(P f )-E nl 



(2-23) 



where Pi(x) is the Legendre polynomial. 

Using this potential, a simple closed-form expression for the S-matrix can be 
obtained, if only one confined and one free channel are considered (see Appendices 
A.1-A.5 of RefP^ for a detailed derivation), viz. 



2a\ 2 Z 



SUE) = l-2i- 



o E{P f ) - E nl 



»fP f j lf (Pfa)h}(P f a) 



(2-24) 



1 + 2iaX 2 £ 



oE(P f )-E nl 



H P f j lf {P f a)hUP f a) 



An exact solution for the 5-matrix can be derived in the most general multichannel 
case as well, resulting in a matrix expression with a similar structure.^ 

The full scattering amplitude can be depicted diagrammatically as shown in 
Fig. HI where the shaded boxes represent the effective meson-meson potential de- 
picted in Fig. [TJ 



MM->MM 



M 



M 



M 



M 



■W777A- 




Fig. 4. The full scattering amplitude for two mesons, with the shaded boxes given by Fig. \T\ 

In order to find resonances in non-exotic meson-meson systems, we search for 
zeros in the denominator of the 5-matrix (Eq. (|2-24p ). which correspond to poles 
in the complex-energy plane. Then we vary the parameters A and a to fit the 
experimental data for a particular case. We also study the pole movements in the 
complex plane, as we change the parameter A, in order to get more physical insight 
into the model, and to evaluate the role of meson loops coupled to quark channels 
in the generation of resonances. 



§3. Perturbative formalism 



The main purpose of this paper is to show that predictions of resonances poles 
based on perturbative calculations in standard quark models can be very misleading. 
In order to demonstrate this in a quantitative way, we now construct a perturbative 
scheme for the formalism described in the previous section. 
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As explained above, the A in our formalism is the coupling of the meson-meson 
•(->• qq vertex. The term corresponding to A 2 thus represents the lowest-order meson- 
meson interaction (meson- meson — >■ qq — >■ meson- meson ). The position of a pole 
in a particular case of meson-meson scattering, above threshold, can be expanded 
perturbatively in terms of A 2 as 



ppole TTi I \2 tpLO , \4 irnNLO ,\6p 



6 171NNLO 



+ 



(3-1) 



The first term of this series corresponds to the confinement pole, which is what we 
should get from the model if the coupling of the quark pair to the meson channels 
vanished. The second term is the leading-order term in A 2 . The denominator of the 
S-matrix (Eq. (j2-24|) ) . written up to leading order in A 2 around the mth confinement 
pole, becomes 



1 + 2iaX' 



9m 



E m + X 2 E m ° _ E 



- fj,kji(ka)hi (ka) = 0, 



which then yields 

E m ° = -^iag m fikji(ka)h}(ka). 
So the pole position (Eq. (|3Tp ). in lowest-order approximation, is 

£P° le f*E m - 2iag 2 m fikj l (ka)h}(ka). 



(3-2) 
(3-3) 
(3-4) 



We now expand the denominator of Eq. (|2-24p to higher order in A 2 . In order 
to do so, we define 

/(A 2 ) = 2iankji(ka)hj(ka), (3-5) 
such that Eq. (|2-24p becomes 



2A 2 /(A 2 ) £ 



9 2 n 



Si(E) 



n=0 



E — En 



On 



(3-6) 



1 + A 2 /(A 2 )E 

n=0 & — 



Now we expand the function / around the energy E = E m , i.e., around A = 0, as 



2 df 

f(E m ) = f(E m ) + A 



d 2 f 



+ 1a 4 

a=o 2 3(A 2 ) 2 a=o ' 6" 9(A 2 ) 3 



a 3 / 



A=0 



+ ... 



f(E m ) + A 



2 df BE 



dE dX 2 



d 2 f ( dE_\ 2 
dE 2 \dX 2 J 



+ 



+ 



df d 2 E 
dE~dX 2 

df d 3 E 
dEdi^f 



A=0 



d 3 f fdE\ 3 d 2 f dE d 2 E 
1 ' +3- • 



dE 3 \dX 2 



dE 2 dx 2 d(x 2 ) 2 



+ 



(3-7) 
(3-8) 



A=0 
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From Eq. (|3-1|) . we have -Em° lc U=o = E m , and 

a ppole 02 P pole 03 ppole 

(J-^m — LO " m „„MT.n f J-J-m 



ox 2 



A=0 



Using these relations in Eq. 



A=o 
we get 

2 1-.LO ®f 



2E„ 



NLO 



9(A 2 ) S 



A=0 



6E; 



NNLO 



f(E m ) = f(E m ) + X 2 E m u ^ 

2 d 2 f 
dE 2 

d 3 f 



E—Er, 



(E L m °) 



+ 6 A 



L0^2 



(-Em ) 



, n pNNLO ^/ 



dE 3 



E — Er, 



E—E„ 

E=E n 
+ ■■■ 



+ 2E„ 



NLO 



dE 



E — En 



1 fi pNLO F LO 9 2 f 
UI2j m ^m QJ£2 



E—Er, 



(3-9) 



It remains to consider the remaining part of the denominator, namely 



n=0 



E — En 



^ 2 E 



<A 



n^rn 



E — Er 



+ A 



2 9n 



E — E m ' 



(3-10) 



which we expand in a series in A 2 , around A = (so at E = E m ), in a similar way 
as the function /, to obtain 



OO 9 

9n 



x 2 y - a 4 £lo y 

4- E m -E n m 4- {E s 



n=0 



E — E n 



9n 



E n ) 2 



m ~ E n f 



.9, 



(E m — E n ) 2 



+9m 



1 



pNLO r^f?NLO\2 pNNLO ^1 



E„f (£LO)2 + \ (£LO)3 ( £LO )2 / 

( , /zr>NLO\3 pNLO ^ pNNLO pN 3 LO 



(E m °y 



(e™Y 



+ ■■■ 



(3-11) 



Upon multiplying Eqs. (|3-9p and (|3-lip . we get an expansion in A 2 for the pole 
position, viz. 



= l + g 2 m E m °f(E,, 



+ A 2 



E 9n (Elf) 
— n 9 m E; 



E — E 
*- J m J -'n 



2 pNLO 

m m 



f(E, 



(E m °) 



-I- n 2 —L 
2 ' 9m a rn 



dE 



E — Er, 



E 

n^m 



9n \E m ) 
(E m — E n ) 2 



+ 9r, 



{(El h °) 



2 pLOpNNLOl I f(E m ) 
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9nE m ° I 9f 



1 



E—E n 



+ 7i9m E ; 



,2 etLO 




E ffn (-^m°) _ \ ffn {E m °) 
(E - E ) 3 ^ (E 
, \ Lj m ■L-'n) \ J - J r, 



E n ) 2 



■2 / / tttNLO\^ opNLOpNNLOpLO , ittN 3 LO /pLO\^\^ fi^i 
9 m \\E m ) ~ E m E m +E m [E m ) j- j -— 



| ^ 9l (^°) 2 , ^ 9 2 n E* LO \ df 

+ I Z> _ R^2 + Z> 



! m En) 2 E m E n 



+ 



E9n { E ln y 2 jpNLO d 2 f 
E m -E n ^ 9m m | dE* 



2 1 ^ K 



<9£ 



E—Er, 



E—E„ 



+ 



1 2 / lq\2 d 3 f 
+ v9m \E m 



6' 



d£ 3 



E—E„ 



Solving this equation order by order in A 2 , we obtain the expressions for e.g. the pole 
position to lowest order, next-to-lowest order, and next-to-next-to-lowest order as 

(3-13) 



E 



m -9mf(Em), 
NLO _ ... I 



= 9 4 m f(E m )%\ E=Em + 9 2 m f(E m ) ^ 91 



E — E ' 



E 



NNLO 



9mf 3 (E m ) ^ Tp _" F n 2 ~~ 9mf 3 (E m ) 



n^m 



E 



9n 



E — E 



39m f 2 (Em) 



' 2^mf {Em) 



E—E„ 



9mf{Em) 



E—E„ 

\dE 



E 



9n 



E — E 



E—E„ 



(3-14) 



Similarly, one can obtain the expressions for even higher-order contributions. 



§4. Results and discussion 



To test the validity of the perturbative calculus, we now choose a few concrete 
examples of meson-meson systems, namely Kir in P wave, DD in P wave, and DK 
in S wave. First we compute the scattering poles in these systems by using the exact 
formalism, explained in Sect. [2j In particular, we study the S'-matrix (Eq. (|2-24p ) 
pole positions in the complex-energy plane, as a function of the coupling A. The 
resulting pole trajectories behave as expected from general considerations: 
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• The unperturbed, or bare, quark-antiquark spectrum is given by Eq. (|2Tip , 
For small coupling to the scattering sector, we must find resonance poles close 
to the levels of this spectrum. 

• For larger coupling, the resonances are expected to acquire larger widths, so the 
pole positions must get larger (negative) imaginary parts. However, associated 
with larger widths are, in general, larger mass shifts. As a consequence, also 
the real parts of the pole positions will deviate substantially from the levels of 
the bare qq spectrum. 

• Sometimes, a further increased coupling may lead to a sufficiently large mass 
shift so as to push the pole below the scattering threshold. For such a situation, 
we expect a bound state. Since a bound state has zero width, the pole should 
eventually end up on the real-energy axis. 

Next, we study the same cases, but now using the perturbative formalism derived 
in Sect. El up to fourth order in A 2 . This way we can compare exact and perturba- 
tive results for the real and imaginary mass shifts due to meson loops. The latter 
quantities are related to predictions for the central resonance masses and resonance 
widths. We will find that, for moderate to large couplings, the perturbative results 
strongly deviate from and do not converge towards the exact ones. In particular, the 
expected behavior of poles below the scattering threshold, namely to move towards 
or along the real-energy axis, cannot be reproduced at all by the studied perturbative 
expansions. 

In all cases we will employ the quark masses m n = m u = rrid = 0.406 GeV, 
m s = 0.508 GeV, m c = 1.562 GeV, as well as the universal oscillator frequency 
ijj = 0.19 GeV, all determined in the model of RefP^ The values of the parameters 
a and A are given separately for each case in the following discussion. 

The parameter a in Eq. (|2-14p describes the average distance at which quark-pair 
creation or annihilation takes place, leading to the two-meson decay of a meson. For 
the case of Kir, which is a nonstrange-strange flavor combination, we choose here 
a = 2.534 GeV -1 . This is smaller than the value a = 3.2 GeV -1 , which would be 
used in a multichannel fit to several mesons. The reason for this discrepancy is that 
one meson- meson channel, namely Kir, has to mimic the effect of a multichannel 
treatment. In order to nonetheless obtain the K*(892) pole at a reasonable position, 
we have to adjust the parameter a. In the DD and DK cases, we use here values 
close to the ones used in a multichannel calculation P^J namely a = 1.72 GeV -1 and 
a = 2.5 GeV -1 , respectively. The price we pay is that, in the DD case, the ip(2S) 
pole does not come out at 3.686 GeV, but about 40-50 MeV higher. 

For the parameter A, which is the universal overall three-meson coupling con- 
stant, we could have used, after scaling, the same value in all three cases. However, 
scaling was not carried ou1^ in the DK case, which makes that the pole now greatly 
overshoots the D s o(2317) position for A = 1. With the correct scaling, it would end 
up roughly 20 MeV too high. 

One might argue that only one value of the coupling A describes the physical 
situation, so that other values are not relevant. However, as analyticity has proven in 
the past to be a powerful tool for constructing scattering amplitudes, the trajectories 
of their poles are also a strong indication for the correctness of their dependence on 
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other parameters. In the following, we will show that perturbative expansions, even 
to higher orders, only have a very limited range of validity, and do not cover the 
realistic case of large couplings^ ,G2DJS2} [ n strong interactions. 

4.1. P-wave Ktt scattering 

The K*(892) resonance is well described by a Breit-Wigner resonance in P-wave 
Ktt scattering, with central mass and resonance width of about 892 MeV and 50 
MeV, respectively. Hence, we expect a pole in the S matrix of Eq. (|2-24p for an 
S-wave nonstrange-strange quark-antiquark system, coupled to a P-wave kaon-pion 
meson-meson system. For the couplings g n we nncP^J in this case 

9 „ = 2 -»(^±l)"\ (4 ,, 

Scattering poles are obtained by studying the zeros of the denominator in the ex- 
pression of Eq. (|2-24p . In Fig. 03(a), we depict the 5-matrix pole positions for a range 
of A values varying from to just over 2. 

In the limit of vanishing coupling, one expects to find the poles at the bare masses 
of the quark-antiquark system, as given by Eq. (|2- . We obtain from Eq. (|2-lip the 
value Poo = 1-199 GeV for the ground-state bare mass, which indeed corresponds to 
the limit of vanishing A along the dashed curve in Fig. [5](a). For larger couplings, we 
observe that the imaginary part of the pole position vanishes at the Ktt threshold. 
This was to be expected, since a large coupling results in a bound state below the 
Ktt threshold, which of course has a zero width. The shape of the pole trajectory 
near the Ktt threshold is in accordance with theory for poles in P-wave scattering 
and also in higher partial waves.^'^ For S-wave scattering the pole behavior is 
different, as we will see in Sect. 14. 3| but again in agreement with theoryP®'^ 

The value A = 1 corresponds to the physical pole, as it roughly reproduces the 
characteristics of the K* (892) resonance. In the present simplified model, the pole 
comes out at (0.972 — «0.026) GeV, as shown in Fig. 03(a). 

The coefficients of the perturbative expansion (Eq. (|3-ip ) are collected in Table HI 
for the case of P-wave Ktt scattering. In Figs. [5](b-f ) we depict the perturbative pole 



coefficient 


value (GeV) 


Eo 

F N°LO 

cJN^LO 
Eo 

Eo 

F N 3 LO 
E 

Eo 


(1.199, 0.) 
(-0.249080686,-0.0878366188) 
(-0.0435913117,0.0471697828) 
(0.0631440181,0.0973648258) 
(0.0869057944,-0.0785515047) 
(-0.067527691 ,-0.0632886834) 



Table I. Coefficients of the perturbative expansion given in Eq. (|3Tp . concerning the pole positions 
of the ground-state pole in P-wave Ktt scattering. 

trajectories for the bare nonstrange-strange qq state at 1.199 GeV. Shown are the 



Drawbacks of applying perturbative schemes to meson spectroscopy 



15 




Fig. 5. The P-wave Kir resonance pole positions in the complex-energy plane, corresponding 
to the second Riemann sheet, under variation of the coupling A. The dashed curves in 
all the plots for N™LO (n = 0, . . . , 4) are the same as the ones shown for N°°LO (labeled 
(a)). The solid curves are the results obtained from the perturbative approximations, 
viz. (b) leading (N°LO, Born) term, (c) next-to-leading (N 1 LO), (d) next-to-next-to- 
leading (N 2 LO), (e) (ncxt-to) 3 -leading (N 3 LO) and (f) (ncxt-to) 4 -lcading (N 4 LO) orders, 
respectively. 



curves for the lowest-order (Born) term (Eq ) and for the next few higher-order 
terms (E^ LO , . . . ), respectively, up to fourth order in A 2 . We find that the Born 
term gives satisfactory pole positions for overall couplings up to A » 0.3. At each 
higher order, the perturbative pole positions, i.e., the central masses and the widths 
of the if* (892) resonance, are better and better determined, up to A ~ 0.75 for the 
fourth-order approximation. However, thereabove things go terribly wrong, and all 
approximations completely fail to reproduce the physical pole at A = 1. 

So we are forced to conclude that perturbation theory is unreliable to describe 
the if* (892) resonance. Moreover, we should add that these higher-order perturba- 
tive calculations are much more tedious than just finding the exact solution for the 
coupled quark-antiquark and meson-meson system. 
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4.2. P-wave DD scattering 

Let us next consider the DD system, which has been studied already a long 
time agopS 1 using the model described in Sect. In RefP^ it was shown that the 
P-wave DD channel, together with higher open-charm channels, can transform the 
bare vector charmonium spectrum into the physical one. In particular, the pole 
stemming from the first radial excitation comes out very close to the '0(25') (3686) 
state, which turns out to contain a significant DD component, besides cc of course. 

The couplings g n in this case are again given by the vector pseudoscalar- 
pseudoscalar vertex, for which we use the same expression as in Eq. (|4-ip . The 
parameter a is now taken at 0.34 fm. Using these inputs, the 5-matrix (Eq. (|2-24p ) 
is calculated, and we search for poles on the second Riemann sheet. We present 
the results of our calculation in Fig. El^a), which depicts the complex-energy plane 
around the mass of the ip{2S) (3686). The dashed line in Fig. EJa) corresponds to 




Fig. 6. P-wave DD pole positions in the complex-energy plane (corresponding to the first 
and second Riemann sheets), under variation of the coupling A. The dashed curves in 
all the plots for N™LO (n = 0, . . . , 4) are the same as the one for N°°LO (labeled (a)). 
The solid curves represent the perturbative results: (b) leading (N°LO, Born) term, (c) 
next-to-leading (f^LO), (d) next-to-next-to-leading (N 2 LO), (e) (next-to) 3 (N 3 LO), and 
(f) (next-to) 4 -leading (N 4 LO) orders, respectively. 



the movement of the pole in the complex plane between the DD threshold and the 
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first radial excitation of the J/0, as A is varied between the limiting values and 1. 
As expected, when the coupling is very small, we find a pole close to the first radial 
excitation of the confined cc spectrum of the model, i.e., near E\ = w(2+3/2)+2m c = 
3.789 GeV (Eq. (|2Tip ). As A is increased to 0.1, the real part of the pole becomes 
Ke(-E') ~ 3.77 GeV Finally, for A ~ 1, the pole is found below the DD threshold, 
very close to 3.7 GeV, which should corresponds to the physical ip(2S)(3G86). 

Next we show the pole positions obtained in the perturbative expansion, viz. 
from Eqs. (|3T2rl3^14")) ) and similar expressions up to fourth order in A 2 . The coeffi- 
cients of Eq. (|3-ip are collected in Table HH 



coefficient 


vaiue (GeV) 


Eo 

i-.N 1 LO 
fro 

fro 

P N 3 LO 
-fro 

pN 4 LO 
-fro 


(3.789, 0.) 
(-0.291265968,-0.0139626856) 
(0.587237403,0.157827196) 
(-0.763203829,-0.81585325) 
(-0.589315239,2.48140259) 
(7.66667227,0.981322456) 



Tabie 11. Coefficients of the perturbative expansion (Eq. (I3-1[) ~1 for the first radiaily excited pole in 
P-wave DD scattering. 

Figure E^b) shows that the pole position found in the leading-order approxi- 
mation agrees with the full calculation only for very small values of A, but as the 
coupling increases, the approximate pole starts to deviate strongly from the exact 
one shown in Fig.[6ja). For example, at A = 0.5, the first-order pole comes out below 
threshold but with a large imaginary part, which is obviously unphysical. For the 
higher-order approximations, the results are even worse, with the pole moving into 
the upper half plane, or extremely deep down in the lower half for the N 2 LO case. 
It becomes evident that no perturbative approximation will produce anything even 
resembling a bound-state pole for A ~ 1. 

4.3. DK S-wave scattering 

Finally, we study the case of S-wave DK scattering taking 

g n = 2- n (n + lf 2 . (4-2) 

Now, as one can see in Fig. [TJ^a), the shape of the pole trajectory near the DK thresh- 
old is very different from the two P-wave cases. For increasing A, the pole approaches 
the real-energy axis below threshold, moves along the axis towards threshold as a 
virtual bound state, and then becomes a bound state, moving finally to lower and 
lower energies. There is a one-to-one relation of this complex-energy pole trajectory 
to the equivalent one in the complex-momentum plane. Thus, a virtual bound state 
moving towards threshold corresponds to a momentum pole moving upwards along 
the negative imaginary axis, passing through the origin when the virtual bound state 
becomes a true bound state, at threshold. In the present case, the pole is still on 
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Fig. 7. iS-wave DK pole positions in the complex-energy plane (corresponding to the first 
and second Riemann sheets), under variation of the coupling A. The dashed curves in 
all the plots for N™LO (n = 0, . . . , 4) are the same as the one for N°°LO (labeled (a)). 
The solid curves represent the perturbative results: (b) leading (N°LO, Born) term, (c) 
next-to-leading (N^O), (d) next-to-next-to-leading (N 2 LO), (e) (next-to) 3 (N 3 LO), and 
(f) (next-to) 4 -leading (N 4 LO) orders, respectively. 



the negative imaginary axis for A = 0.4, but already on the positive one for A = 0.5. 
This phenomenon, which happens exclusively for S-wave scattering, as can be seen 
from the effective-range expansion, is well described in RefsP^ For A » 0.6, the 
bound-state pole reproduces the -D<jo(2317) mass. 

The coefficients of Eq. (|3Tp for the case of 5-wave DK scattering are collected 
in Table [TTTl One sees at a glance that these coefficients do not promise any kind of 
convergence. Indeed, upon inspecting Figs. [7|(b-f), one notices that only for small 
values of A the perturbative pole positions agree with the exact ones. However, for 
A > 0.3, the discrepancies grow rapidly, and no significant improvement is observed 
for higher orders of perturbation theory. 



§5. Summary and conclusions 



We have studied the discrepancies between perturbative estimates for resonance 
pole positions and the exact ones, in the context of a simple soluble model for 
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coefficient 


value (GeV) 


Eo 
Eo 

zJ^LO 
E 

Eo 

Eo 

Eo 


(2.545,0.) 
(-0.445872986,-0.67333385) 
(-1.36316635,-1.84200144) 
(-9.95402765,-8.57593239) 
(-68.4326606,-43.6873369) 
(-299.284654,-138.032657) 



Table III. Coefficients of the perturbative expansion (Eq. (OS}) for the ground-state pole in S-wave 
DK scattering. 

hadronic decay of a meson. In none of the considered cases satisfactory results were 
obtained with the perturbative method, and not even any significant improvement 
was found for increasing orders of perturbation theory. In particular, for bound states 
below the lowest strong-decay threshold, no perturbative approximation produced 
anything like a pole close to the real-energy axis. But also in the case of a normal 
and not even very broad resonance, namely the if* (892), the perturbative approach 
failed completely. 

These results should be a warning for quark-model builders, because of two 
reasons. First of all, the found large real mass shifts are, as a consequence of an- 
alyticity, inseparably connected to the generation of the physical hadronic widths, 
as demonstrated here for the if* (892), but shown already many years ago for a va- 
riety of mesons,^ and confirmed in several more recent papers referred to above. 
Therefore, any spectroscopic conclusions based on single-channel, "quenched" quark 
models should be taken with a great deal of caution. The second reason is that even 
those quark models which pay some attention to strong decay, usually do this by 
employing perturbative methods. The results presented here make it clear that a 
completely non-perturbative treatment of hadronic resonances and bound states is 
required for a realistic description. 
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